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1 INTRODUCTION 



The intergalactic medium (IGM) con tains most of th e 
cosmic baryons at all epochs (e.g., iDave et al.l l200l[ ). 
It also 
als (e.g., 



|w iersma et al.. .2009b. ) , products of star formation within 
galaxies that have been dispersed out o f galaxies presumably 
via large-scale galactic outflows (e.g.. lAguirre et al.ll2001al : 
lOppc nhcimcr & Dave 2006). The distribution and physical 
state of these metals thus provide a powerful tracer of galac- 
tic outflow processes, which are a key uncertainty in current 
galaxy formation models. 

[Oppenheimer et all l|2012D argued from cosmological 
simulations with outflows that the IGM is enriched in an 
"outside-in" manner, with low-density regions enriched at 
earlier epochs, and that subsequently produced metals tend 
to remain closer to gala xies over the last 10 G y r of c osmic 
time. In simulations by lOppenheimer fc Pavel l|2008h that 



ABSTRACT 

We study the physical conditions of the circum-galactic medium (CGM) around 
z = 0.25 galaxies as traced by Hi and metal line absorption, using cosmological 
hydrodynamic simulations that include galactic outflows. Using lines of sight targeted 
at impact parameters from 10 kpc to 1 Mpc around galaxies with halo masses from 
10^^ — IO^^Mq, we study the physical conditions and their variation with impact pa- 
rameter b and line-of-sight velocity Aw in the CGM as traced by Hi, Mgii, Siiv, 
Civ, Ovi, and Neviii absorbers. All ions show a strong excess of absorption near 
galaxies compared to random lines of sight. The excess continues beyond 1 Mpc, re- 
flecting the correlation of metal absorption with large-scale structure. Absorption is 
particularly enhanced within about A?; < 300 kms~^ and roughly 300 kpc of galaxies 
(with distances somewhat larger for the highest ion), approximately delineating the 
CGM; this range contains the majority of global metal absorption. Low ions like Mgii 
and Si iv predominantly arise in denser gas closer to galaxies and drop more rapidly 
with &, while high ions Ovi and Neviii trace more diffusely distributed gas with a 
comparatively flat radial profile; C IV is intermediate. All ions predominantly trace 
T - IQ-^-^'-^K photo-ionised gas at all 6, but when hot CGM gas is present (mostly 
in larger halos), we see strong collisionally- ionised Ovi and Neviii at 5 ^ 100 kpc. 
Larger halo masses generally produce more absorption, though overall the trends are 
not as strong as that with impact parameter. These findings arise using our favoured 
outflow scalings as expected for momentum-driven winds; with no winds, the CGM 
gas remains mostly unenriched, while our outflow model with a constant velocity and 
mass loading factor produce hotter, more widely dispersed metals. 

broadly reproduce IGM enrichment observations, outflow- 
ing metals reach a median distance of roughly 100 phys- 
ical kpc before turning around, which is well outside the 
virial radius of small early galaxies but results in "halo foun- 
tains" within larger, present-day systems. Hence by low red- 
shifts, a large fraction of the ejected metals are expected 
to reside close to galaxies, within what is now frequently 
called the circum-galactic medium (CGM). Indeed, obser- 
vations indicate that C iv, Mg ii and O vi absorbers can 
often be asso ciated with galaxies with i n tens or hundreds 
of kpc (e.g., ISteidel fc Sargent! Il99j: IStocke et all l2006l: 



c ontains a substant i al fraction of cosmic met- 
Ferrara ct al.j l2005l : iDave fc Oppenheimer! 1200?! : 



IWa kker fc Savagel 120 091: 'Pro chaska et aLlbOllT I Chen et al! 
200l': IChen fc Mulc hacv 200!^). The CGM thus provides a 



fossil record of early IGM enrichment combined with the 
gravitational growth of structure, in addition to metals de- 
posited by more recent outflows. 

Recent observations using Nubble's Cosmic Origins 
Spectrograph (COS) provide a g reatly enhanced vi ew o f 
CGM metals at low redshifts. iTumlinson et al! (|201ll ) 
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showed that the CGM around star-forming galaxies con- 
tains an amount of oxygen comparable to or possibly exceed- 
ing that in the interstellar medium (ISM) of those galax- 
ies, thereby directly demonstratin g that metals have been 
expelled en masse from the ISM l|Tripp et alj boill ). The 
sensitivity of COS not only permits the selection of more 
numerous, fainter background sources that can probe the 
CGM of low-redshift galaxies, but also enables the detec- 
tion of weaker and less common ions such as Nevill and 
Si IV. These data will provide new and stringent constraints 
on theoretical models of outflows that connect galaxies and 
their surrounding gas. 

Cosmological hydrodynamic simulations offer a 
uniquely well-suited platform to study the co-evolution 
of galaxies and their surrounding gas. Such simulations 
dynamically incorporate galactic outflows driven by star 
formation, thereby self-consistently enriching the IGM and 
CGM with the by-products of stellar death, as well as 
accounting for the movement and physical state of such 
metals driven by hierarchical structure formation. With 
simulations of representative cosmological volumes, one can 
test outflow models by comparing the predicted statistics 
directly to absorption line observations of Hi and metals, 
around galaxies and in the diffuse IGM. Tracking the 
complex inflow and outflow processes that govern both 
galaxies and the IGM is necessary to fully understand how 
IGM and CGM enrichment occur, and to properly interpret 
present and upcoming data with COS and other facilities. 

In this paper we continue our exploration of Hi and 
metal absorption in the low-redshift IGM using cosmolog- 
ical hydr odynamic simulat ions, building on the work pre - 
sented bv lDave et al] l|2010l ') and lOppenheimer et al] (|2012h . 
In those papers we studied absorption properties along ran- 
dom lines of sight through cosmological volumes, which only 
occasionally intercepted the CGM of galaxies. In this work, 
we select targeted lines of sight at various impact parameters 
around galaxies in our simulations to directly probe the ab- 
sorption properties and physical conditions of gas within the 
CGM. The goal of this paper is to understand the relation- 
ship between galaxy properties and HI and metal absorbers 
that probe the gas around those galaxies. 

We explore a range of metal ionisation potential states 
from low (Mgll) to high (Nevill), as well as Hi, and find 
that the behaviour of many CGM absorption line proper- 
ties shows distinct trends with the ionisation level of the 
tracer species. CGM absorption lines also show trends with 
halo mass. These predictions set the stage for a more de- 
tailed understanding of CGM observations and provide a 
more complete picture for how metals trace the motion of 
gas in and out of galaxies. 

Our paper is organised as follows: In !j2]we explain our 
methods, in |3]we show the density and temperature of ab- 
sorbers, in 52] we explore absorption around galaxies in red- 
shift space, in !j5]we study absorption as a function of impact 
parameter, in |6] we examine the column density distribu- 
tions as a function of impact parameter and halo mass, in JT] 
we examine trends with outflow model, and in [J8]we present 
our conclusions. 



2 SIMULATIONS & ANALYSIS 



2.1 The Code and Input Physics 



We use our modified version (|Oppenheimer fc Day j [20081 ) 
of the N-bo dy-|-smooth pa rticle hydrodynamic (SPH) code 



Gadget- 2 (ISpringell 



in §2.1 of iDave et al 



20051'). which is more fully described 
(120101 '). The only code update since 



lOppenheimer fc Dave (|2008|) is the option to include metal- 



line cooling rates from lWiersma et alj||2009al ). as we discuss 
below. 

Briefly, Gadget-2 computes gravitational forces on 
a set of particles using a tree-particle-mesh algorithm 
and uses an entropy-con serving formulation of SPH 
jSpringel fc HernauistH2002l ') to simulate pressure forces and 
hydrodynamic shocks. We include radiative cooling assum- 
ing ionisation equilib rium for primordial species foUow- 
ing Katz et al.| (119961) and metals based on the tables of 
|w iersma et al.l (l2009a|) that assume ion isation equilibrium 
in the presence of the lHaardt fc Madaul (|200ll ') background. 
In our older simulations, and our ancillary wind model sim- 
ulations in this paper, we employ metal-line cooling rates 
based on the coUisi onal ion i sation equilibrium (CIE) mod- 
els of ^Sutherland fc Dopital (| 19931 '). The latter rates incor- 
rectly over estimate metal-line cooling as photo-ionisation 
equilibrium (PIE) s uppresses coo ling of metal-enriched gas 
as dem onstrated by ISmithI (|201ll ) and iTepper-Garcfa et al] 
(|201ll ). We show comparisons of these two cooling models 
iniZl 

Star formation follows a ISchmidtl (| 19591 ) Law cal- 
ibrated to the Kennicut3 (| 19981 ') relation, following 
ISprineel fc HernauistI (|2003l ). The interstellar medium 
(ISM) is modelled using the analytic sub grid recipe of 
iMcKee fc Ostrikeil (| 19771 ). where supernova energy is re- 
turned to ISM particles using t he tw o-phase SPH formu- 
lation of ISpringel fc Hernauisd (|2003l ). Star particles are 
spawned from ISM particles probabilistically according to 
the instantaneously calculated star formation rate (SFR); 
an ISM particle can spawn up to two star particles. 

Star formation-driven kinetic feedback is implemented 
in these simulations by giving a velocity kick {v^) to ISM 
particles chosen probabilistically at a rate proportional to 
their star formation rate. The ratio of the mass outflow rate 
to the star formation rate is termed the mass loading factor 
{rf). We use the relat ions for v-„ an d t? of momentum-driven 
winds formulated bv lMurrav et al.l (|2005l ): 



= 3cr%/fr 
= CTo/cr, 



(1) 
(2) 



where f^ is the luminosity in units of the Eddington lumi- 
nosity required to expel gas from a galaxy potential. Go ~ 
150 kms~^, and a is the gala xy's internal velocity dispersion 
(|Oppenheimer fc Davel200^ ). We refer to this wind model as 
"vzw" and refer the reader to lOppenheimer fc Pavel (|2008l ) 
for a more complete description. We also evolve a model 
without winds ("nw") and a constant wind ("cw") model 
with iiw = 680 kms~^ and r) — 2 for all galaxies, as de- 
scribed in iDave et all (|2010l) . 
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Table 1. Comparison of r48n384vzw models. 



was feat ured prominently in a number of our previous publi- 
cations ( Oppenheimer et al .I2OIOI : iDave et al...2010. . ,2011a ibl: 



Z-Cooling 


CIE° 


PIE'' 


/bar(T ^ 105K)'= 


0.394 


0.408 


/metals(T ^ lO^K)'^ 


0.056 


0.066 




0.084 


0.076 


PSFR. 


0.074 


0.066 


dN/dz(7Vjj, > 10"cm-2)9 


8.88 


7.94 


oft 


57.6 


38.8 


nh 
»'CIV 


33.9 


23.7 


Qh 

"SilV 


10.2 


7.9 



"CoUisional ionisation equil ibrium metal cooling rates from 
[Sutherland fc Dopital l ll993l'l . 

''Photo-ionisation equi librium metal cooling rates using 
IWiersma et al.l l|2009al '). 

'^Fraction of diffuse, non-ISM baryons with T > 10^ K at 2 = 0.25. 
''Fraction of diffuse, non-ISM metals with T > 10^ K at z = 0.25. 
"^Fraction of baryons in stars at 2: = 0.25. 
^Instantaneous SFR at 2: = 0.25 in Mq yr~^ Mpc"^ 
^Frequency of Hi components with A^jj j ^ 10^*cm~^ at 2 = 0.25. 
''f2 summed from all gas outside galaxies a,t z = 0.25, in units 
of 10-8. 



2.2 Simulation Runs 

Our main model in this paper is a 48/i~^Mpc, 2 x 384'^ sim- 
ulation with vzw winds. This simulation is identical to the 



r48n3 84vz w simulation used by Oppenheimer et al.l (|2O10l . 
20121) and iDaye et al.) (|201(]| . l2011al lb^ except that it uses 



Wiersma etaP (|2009al ;) metal-line cooling rates. 



We choose the vzw wind model as our fiducial model 
owing to its successes fitting important properties of galax- 
ies and the IGM at a variety of redshifts. At high red- 
shift, simulations with vzw winds provide adequate fits 
to th e z = 2 ma s s-meta llicity relationship (jFinlator et al.l 

l2008l : iDave et al.l l2011a|). the z ^ 6 galaxy luminosity 

function ( Dave fc Oppenheimeill2007l : iFinlator et a l 203JL) , 
and obs ervations of IGM metal enr i chment at redshifts 

> 1.5 (|Oppenheimer fc Pavel l2006l : lOppenheimer et al] 

|2009| '1. Simulations evolved to 2 = r eproduce the ob- 
serve d gala ctic stellar mass funct ion (Sell et al.l l2003l : 



Baldrv. Glazebrook. fc Driver 1 1200^ ) below 5 x lO" M(, 



I Oppenheimer et al. 20101). variou s statis tical properties of 
present-day galaxies (|Dave et al.l l2011al lbh. and Hi and 
met al-line statistics from quasar a.bsorption line (QAL) spec - 
tra llOppenheimer fc Dayg|2009l : lOppenheimer et alll2012h . 
lOppenheimer et al.l (|2012l ) explored a small, 16h ^Mpc, 2 x 
128^ sim ulation incorporating vzw winds and lWiersma et aU 
l|2009al ) metal-line cooling, finding it to be their preferred 
model owing to its theoretically motivated feedback pre- 
scription capable of reproducing the above-listed observa- 
tions, co mbined with a more correc t treatment of metal-line 
cooling. lOppenheimer et al.l (|2012l ') also compared the vzw 
model, as well as a constant wind and no wind model, to 
observations of random lines of sight, as shown in Figures 
6-8, 10, 14-15, 18 and 21. We note that these simulations do 
not include any additional mechanisms such as black hole 
feedback to quench star formation in massive g alaxies to re- 
prod uce the observed galaxy red sequence fe.g.. lGabor et ah! 

Because the old r48n384vzw evolved with CIE cooling 



lOppenheimer et al.ll2012l ). we list the differences in some rel- 
evant quantities relative to using PIE cooling rates in Table 
[T]at z = 0.25. It shows that the differences that result from 
replacing CIE with PIE metal-line cooling rates are minor 
to moderate. The fraction of baryons in diffuse (i.e. non-star 
forming) gas with T ^ 10^ K increases by 3.6%, while the 
associated metals increase by 18%. Less cooling allows more 
baryons and especially metals to remain in the warm-hot 
intergalactic medium (WHIM), defined as where T > 10^ K 
and 5 < 5th, where Sth is the division between bound an d 
unbound gas (lOppenheimer et al.|[2oT^ : iDave etHI 120101 ). 
The total integrated star formation (i.e. baryon fraction in 
stars) is reduced by 10%, and the instantaneous star forma- 
tion rate is reduced by 11%, both because enriched winds 
are less capable of cooling an d hence recycling onto galax- 
ies (|Oppenheimer et al.ll2010l ). Considering Lya absorption 
tracing the IGM, the frequency of single Voigt-profile-fitted 
components with A'jjj ^ lO^^cm"^ drops by 10.6%. Com- 
mon metal ion species observed by COS all decline by 23 - 
33%, which was also found by lOppenheimer et al.l (|2012l ) 
and shown to provide as good if not better agreement with 
the available observations. Hydrogen and metal-line absorp- 
tion primarily arise from T < 10^ K gas, which is reduced 
with PIE cooling, and more so for metals. The relative differ- 
ences between these two vzw simulations using CIE or PIE 
metal cooling is less than it would be for a stron ger wind 
model such as that used in Wiorsma et al.1 l|2009ah . because 
with vzw winds diffuse metals are placed at higher over- 
densities where they cool rapidly in either cooling scheme 
lOppenheimer et aLll2012l ). 

In §7, we compare these vzw wind models to two other 
wind models, a constant wind and a no wind model, both 
with 48/i-^Mpc size and 2 x 384 ^ resolution. The con - 
stant wind model presented here has I Wiersma et al.1 (l2009ar) 
mctal-line cooling, included because Oppenheimer et al.l 
(2012) found that PIE cooling was more important for met- 
als pushed to lower densities by these stronger winds. We 
use th e same no wind simulation as in lOppenheimer et al.l 
(|2010t ). Since metals rarely reach densities where PIE metal 
cooling is important, there is no need to rerun a no wind 
simulations with PIE metal cooling. 

Finally, the cosmology used in our r-series simulations 
is: fim = 0.28, = 0.72, h = 0.7, spectral index n = 0.96, 
o-g = 0.82, an d flb = 0.046. This agrees with the WMAP- 
7 constraints (|jarosik et al.l I2OIII ) . The gas particle mass 
is 3.56 X 10^ Mq, and the dark matter particle mass is 
1.81 X IO^Mq giving an effective galaxy mass resolution of 
about 2.3 X IO^Mq and a dark matter halo mass resolution of 
about 1.21 X IO^^Mq. In this paper we will focus on galaxies 
and halos generally well above these resolution limits. Un- 
less otherwise specified, we quote distances in physical (not 
comoving) units, and we usually refer to them in Mpc rather 
than Mpc. 

2.3 Generating Spectra with Specexbin 

We use our spectral generation code Specexbin to calcu- 
late physical properties of the gas. Specexbin is described in 
more detail in i;2. 5 of,Oppenheimer fc Da ve (2006) and more 
recently in §2.3 of lDave et al.l ( 201ol ). Briefly. Specexbin av- 
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erages physical properties including the gas density, the tem- 
perature, the metallicity and the velocity of SPH particles in 
physical coordinates along a sight line. It then uses loo k-up 
tables calculated with CLOUDY jPerland et al.lll998l . ver- 
sion 08.00) to find the ionisation fraction for the relevant 
ionic species. 

These look-up tables, functions of density, tempera- 
ture, and redshift, have been calculated assuming coUi- 
sional plus photo-ionis ation equilibrium with a uniform 
iHaardt fc Madaul background. Specexbin then con- 

verts to velocity coordinates using Bubble's Law, account- 
ing for the peculiar velocities of the SPH particles, and 
adds thermal broadening using the temperature and atomic 
weights of the various ion species. To match the mean ob- 
served absorption we multiply the Haardt fc M adau (2001,) 
spectrum by factors of 1.5 for the vzw model, and 1.11 for the 
no win d and constant wind models, as found in lDave et al.l 
(|2010l ). 

Since our last published work with Specexbin, we have 
updated the code to include an approximate physically mo- 
tivated self-shielding from the ionisation background. This 
is significant here owing to the fact that we are interested 
in absorption close to galaxies that may arise in fairly dense 
(and possibly self-shielded) gas. The self-shielding correc- 
tion is purely local and is done on a particle-by-particle ba- 
sis. We compute the fraction of each particle's mass that is 
internally self-shielded from the ionising background by in- 
tegrating the Hi column density inwards from the particle's 
edge, assuming that it has a density profile given by the SPH 
smoothing kernel, until the column density crosses a thresh- 
old; within that radius we assume that all the hydrogen is 
neutral. We choose a threshold of 10^*cm~^, which provides 
the best fit to the neutral fractions obtained in the ful l ra- 
diative transfer models of iFaucher-Giguere et al.l (|2009l '). In 
practise, because of a reasonably tight c orrelation betwee n 
column density and physical density (e.g.. lDave et al.ll201(J ). 
this results in a fairly sharp density threshold of 0.01 cm~^ 
above which the gas is fully neutral. If the gas density is 
greater than 0.13 cm~"^ our star formation algorithm as- 
sumes that it is star-forming. In this case we assume that 
the gas is completely neutral, which is reasonable since the 
cold fraction in the two-phase medium always dominates by 
mass over the hot fraction. This choice of self-shielding cor- 
rection also moves all the magnesium into Mgll, which is 
what one expects for magnesium in neutral hydrogen based 
upon the ionisation potentials of magnesium. For random 
sight lines, we rarely probe these densities, but for targeted 
lines of sight particularly within 10 kpc of a galaxy, the ef- 
fects of self-shielding can be significant. We also put in a 
tra nsition from neutral to mo lecular HI for ISM gas using 
the lBlitz fc Rosolowsk\^ l|2006l ) pressure criterion. 

Our spectral extraction implicitly assumes that all met- 
als are in the gas phase, but the observed correlation of 
quasar colours with projected separation from foreground 
galaxies sugges ts sub s tantial amounts of intergalactic dust 
(,Menard et aLlbOld '). IZu et al.l (l201ll) show that o ur vzw 
simulations can reproduce the llVlenard et al. I ()2010l ) obser- 
vations if roughly 25% of extragalactic metals (by mass) are 
depleted onto dust with an SMC-like grain size distribution. 
Precise accounting is difficult because of uncertainties in the 
observations, in the composition and size distribution of the 
dust grains, and in the destruction and production rates as 



a function of time and galactocentric distance. Dust deple- 
tion is an inevitable source of uncertainty in metal-line ab- 
sorption predictions from cosmological simulations; it could 
plausibly lower our predicted absorption for refractory ele- 
ments such as carbon and silicon by 25-50%. 

2.4 Ion Selection 

We investigate six different species: Hll216, Mgll2796, 
Siivl394, Civl548, Ovil032, and Nevill770. These repre- 
sent some of the most common metal ions that COS probes, 
along with Hi. All the metal lines have doublets, making 
their identification in observed spectra more straightfor- 
ward. While Hi in the Lya forest is associated with diffuse 
baryons tracing large-scale structures in the Universe (e.g., 
iDave et al.ll 19991 ) , stronger H I absorbers ar e associated with 
higher-density structures in the CGM fe.g.. lFumagalli et al.l 
I2OIII : I van de Voordlionl ) and are sensitive to the outflow 
model used in the simulation (e.g.. lDave et al.ll2010t l as well 
as the cooling implementation (see Table [l]). 

Mgll is one of the most common observed lines in 
QAL spectra, and it appears to be associated with the 
CGM envi ronments of ga l axies , especially star-forming 
ones (e.g., iKacprzak et al.l 20071), and ma y indicate re- 
cent o utflows fe.g.. | Bordoloi et al.l I2OIII : iBouche et al.l 
l2012l: [Mateiek fc Simcod |2012| ) and/or inflowing gas 
teacprza k et al.ll2010l )~ Mg II is often associat ed with H I ab- 
sorbe rs of A^jjj = 10"-^ - 10^^ °cm-^ (e.g., iHellsten et aD 
Il998l ). indicating that Mgll frequently arises from self- 
shielded gas. Since our cosmological simulations do not have 
the resolution to fully resolve this self-shielding, we make 
the assumption that all magnesium in a self-shielded re- 
gion is Mgll, because the ionisation potential of Mgll is 
above the Hi ionisation energy while Mgl is below it. We 
cannot resolve the true substructure of Mg 11, but we hope 
to understand the typical densities and environments from 
which this absorption arises. This prescription produces 
eq uivalent widths tha t roughly match recent observations 
bv lRubin et al.l (|2012l ). suggesting that our Mgll modelling 
is not dramatically wrong. 

Si IV and C iv are relatively strong lines that fall within 
the far-UV window often probed by COS. These are mid- 
ionisation potential lines that typically arise in optically 
thin g as, but still probe fairly ov erdense regions at low red- 
shifts (jOppenheimer et al.ll2012l ). 

Ovi is the most commonly observed high 
ionisation po tential metal absorb er at low red- 
shifts (e.g., Tumlinson et a l.l I2OIIJ; lOanforth fc Shulil 



iThom fc Cheiil l2008l: Brcgmai^ l2007l: IStocke et al 
'Ch en fc MulchaevI |2009(: Proch aska et aJj 
Sembach et al.l l2003l: iTripp et al l 1200 0, 200^)^ 



Oppenheimer fc Davel l |2009l ) examined the nature of diffuse 
Ovi absorbers in these simulations using random lines 
of sight through their simulation volume and found that 
it traces mostly photo-ionised gas in the diffuse IGM, 
warm-hot gas, and gas within halos. Here we focus on O vi 
absorption within ha lo gas, which we wi ll comp are to the 
previous results from lOppenheimer et al.l l|2012l ) on diffuse 
IGM Ovi. 

The highest ionisation potential line that we consider is 
Nevill, and it is va luable because it can trace 10^"^ K gas 
(jSavage et al.1 12005| ) better than any other UV resonance 
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line. Nonetheless, as argued in lOppenheimer et al.l (|2012h . 
very weak Nevill can arise in gas photo-ionised by the 
metagalacti c back grou nd. Recent observat i ons o f Ne vill by 
iTripp et all (HoTl) and lMulchaev fc ChenI (|2009l ) have been 
interpreted to suggest that Nevill absorption arises from 
a transitional phase on the surfaces of cold clouds moving 
through hotter material. Unfortunately, our simulations are 
not able to resolve such surfaces, so we may be missing this 
component of Ne vill absorption in our models. Nonetheless, 
we do find a significa nt number of such absorbe rs, in rough 
agreement with data (jOppenheimer et al.ll2012l ). 

2.5 Generating and Analysing Simulated Spectra 

We examine both targeted and random lines of sight (LOS). 
For our targeted LOS we randomly select central galaxies 
for each of three difi'erent halo mass bins: ioi° '^5-ii.25 
(labelled 10^^ Mq) ; ioii-75-i2.25 (labelled 10^^ Mq) 
and 1012-75-13.25 (labelled 10^=^ M©). For lO^M© and 
IO^^Mq, we select 250 galaxies, while for lO^'^M© there are 
only 86 central galaxies in the simulation so we use all of 
them in our sample. Because we do not properly model the 
detailed internal structure of the ISM in our galaxies, we do 
not present results through the centres of the galaxies, but 
rather start at 6 = 10 kpc. We choose impact parameters b 
ranging from 10 kpc (centre) out to 1 Mpc, with the spacing 
increasing with b. For each impact parameter b, we produce 
four LOS per galaxy, meaning x+b, x-b, y+b, y-b, for a total 
of 1,000 LOS per b per mass bin (344 for W^^Mq). 

We also generate random LOS, to compare with our 
targeted LOS. For our random LOS, we choose a 100 x 100 
grid spanning the simulation in x — y space, for a total of 
10,000 LOS. These are sim ilar to the spectra generated for 
lOppenheimer et al.l l|2012l '). 

After generating optical depths with Specexbin, we 
construct artificial spectra by convolving our data with 
cos's line spread function (LSF) and adding Gaussian ran- 
dom noise with S/N=30 per 6 kms^^ pixel. The LSF is 
roughly Gaussian with a ~ 17 km/s FWHM, but it has some 
non-Gaussianity in the wings. Also, there is a small wave- 
length dependence to the LSF; for simplicity, we choose the 
G130M LSF at 1450A as representative of all wavelengths 
to avoid complications in our interpretations owing to vari- 
ations in the LSF. 

Figure [T] shows examples of our simulated spectra. For 
each column, a single central galaxy is chosen from that 
mass bin, and we plot spectra at impact parameters of 10 
kpc (red), 100 kpc (green), and 1 Mpc (blue) away from 
that galaxy. These particular galaxies were chosen for this 
figure since they have many absorption features, which illus- 
trates the methodology well, although that is not necessarily 
typical. 

The top two panels show the gas temperature and den- 
sity weighted by Hi optical depth, and the lower panels 
show simulated spectra for Hi and other ions as labelled. 
As expected, the gas density is higher at smaller impact 
parameters, which yields stronger absorption in most ions. 
Also, more massive halos generally have more absorption. 
While this figure only shows three selected galaxies, one can 
already see some of the trends that will be explored and 
quantified later in this paper, although other trends are sta- 
tistical in nature and not well depicted in this figure. 



We fi t Voigt profiles t o the absorption features using 
Auto VP ^^^^\W9^), which yields column densities, 
line widths, (rest) equivalent widths, and redshifts for each 
absorber. We set the detection significance criterion to 4(t, 
which should be quite conservative given that the noise is 
Gaussian random by construction. Since the details of de- 
blending absorption features into multiple components can 
be quite sensitive to the noise level and algorithmic details, 
we generally focus on absorption systems, in which we com- 
bine all lines that have separations < 100 km/s. This mit- 
igates the sensitivity to these issues, and allows for a more 
robust comparison to the observations, at the cost of dis- 
carding some of the information available about the internal 
kinematics of the absorption features. Some of our figures 
below plot median line properties, and since weak lines out- 
number strong lines, median properties necessarily depend 
on the lower cutoff. Our results should be taken to refer 
to the population of lines detectable with typical S/N=30 
COS spectra, with closely separated components combined 
into systems. 



3 PHYSICAL CONDITIONS OF ABSORBERS 
3.1 The CGM in Absorption 

We seek to understand the density and temperature of the 
CGM and how those conditions give rise to the presence of 
various observable ions. We begin by presenting a pictorial 
overview of the simulated CGM, as seen in physical condi- 
tions and line absorption, to build our intuition as to how 
these physical conditions relate to observables. 

Figure [2] shows a stacked median image of the galaxies 
around which we generate our targeted lines of sight. We 
generate the images by taking the pixel-by-pixel median of 
each given quantity for the 250 galaxies in the mass bin of 
10^1 M0 (top row and IO^^Mq (middle row), and the 86 
galaxies in the 10^^ Mq mass bin (bottom row). Because this 
is a pixelized median, the satellites are effectively removed. 
The depth of the image is 658 kpc centred at the redshift 
of the galaxy (~ 1.2% of the full simulation depth). The 
first two columns show the median overdensity and median 
temperature of all the gas. The remaining plots show, in each 
column, the absorption in the various ions that we consider 
in this paper, again for each of our three halo mass bins. 
Note the colour scale for HI, subdividing Hi into damped 
Lymana systems(red) , Lyman-li mits sytems (gre e n) an d the 
Lyg forest (blue), consist ent with lPontzen et al.l l|2008l ) and 
Ivan de Voort et"^ l|2012l ). is different from the metal lines. 
We also include, in the lower two panels, results for two 
different wind models, a constant wind "cw" model and a 
no wind "nw" model for the 10^^ Mq halo mass bin (see §7). 

As expected, the gas density is strongly concentrated 
towards the (stacked) central galaxy. Galaxies in more mas- 
sive halos have a larger extent. The temperatures increase 
with mass bin, with the W^^Mq mass bin showing halo 
gas that is predominantly near the virial temperature, but 
other halos showing more sub-virial temperatures. This re- 
flects the now well-understood division at IO^^Mq between 
cold and h ot-gas-dominated halo s (e.g.. [ Keres et al. I l2005l : 
iDekel fc Birnboim 2006 ; Gabor etal.l 120111 ). Some cool gas 



must be present in 10 Mq mass halos, however, since we 
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Figure 1. Simulated spectra for tliree galaxies at z=0.25, convolved with the COS line spread function and noise added with a S/N=30. 
Each column corresponds to a representative galaxy with the labelled halo mass. We plot lines of sight at distances of 10 kpc (red), 100 
kpc (green), and 1 Mpc (blue) away from the galaxy. The top two panels show the Hi optical depth- weighted density and temperature, 
and the lower panels show simulated spectra for Hi and other ions. All units are physical, for a Hubble parameter h = 0.7. 



find significant Hi absorption (see §3.4). We reiterate that 
our simulations do not include a prescription to quench star 
formation in hot halos, which is required to produce red and 
dead galaxies (Gabor et al. 2011), so more realistic models 
might have less cool gas in IQ^^Mq halos. However, recent 
observations (Thom et al. 2012, in preparation), suggest that 
even early type galaxies have prevalent Hi absorption, so 
perhaps whatever mechanism quenches star formation does 
not greatly impact the cooler halo gas. 



The morphology of the absorption (neglecting individ- 
ual satellite contributions) is strongly dependent on the ion 
being probed. Mgll is highly centrally concentrated, and 
shows very little qualitative morphological difference be- 
tween the different halo masses (on this plot scale). At the 
opposite extreme, Nevill shows very diffuse, extended ab- 
sorption in the lower mass halos, and is more centrally con- 
centrated in the highest halo mass bin. The other ions lie in 
between these two extremes, progressing smoothly through 
Si IV, C IV, and O vi. We foreshadow our upcoming results by 
ordering these ions by increasing ionisation potential, which 
we will demonstrate provides an intuitive guide towards un- 
derstanding how absorption traces physical conditions. 



3.2 Phase Space Plots 

To begin to quantify the physical conditions giving rise to H I 
and metal absorption around galaxies, it is useful to examine 
the location of absorption systems in density-temperature 
phase space. This provides a valuable overview of the phys- 
ical conditions of the gas traced by the various absorption 
species and sets the stage for discussing how various observ- 
able properties trace these underlying physical conditions. 
We thus begin by examining where our systems appear in 
phase space as a function of impact parameter. 

In Figures |3HS1 we show the overdensity-temperature 
space locations of absorption systems (circles, colour-coded 
by column density) at different impact radii as labelled, over- 
laid on the global (i.e. cosmic) mass distribution of that 
ion (fiion) computed for all non star- forming gas {jih < 
0.13 cm~^) in the simulation volume (grey shading). We plot 
this for galaxies in our 10^^ Mq halo mass bin. The verti- 
cal dashed line shows the cosmic mean density at z = 0.25 
(riH ~ 10"® '*'' cm~^), and the solid horizontal and verti- 
cal lines subdivi de this cosmic pha se space into the four re- 
gions defined bv lDave et all (|2010l ). The temperature cut at 
T = 10^ K divides "hot" phases from "cool" phases, and the 
density cut distinguishes halo gas (riH > 10"'* '*^ cm"'', i.e. 
the virial radius density at z = 0.25) from diffuse, intergalac- 
tic gas. The four quadrants, going clockwise starting at the 
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Figure 2. Columns one and two show the median log overdensity and median log temperature; subsequent columns show median log 
column densities for each ion as labelled, all at z=0.25. First three rows are for vzw; the first row is for mass bin IO'^^Mq, second 
for IQi^ Mq, and the third for 10^^ Mq. The fourth and fifth rows are for 10^'^ Mq halos in the constant wind and no wind models, 
respectively, and are discussed further in §7. Note the colour scales for Hi and Nevill are different than for the other ions. All panels are 
658 kpc across. 



upper right, therefore, correspond to: hot halo, condensed, 
diffuse, and WHIM gas. 

In targeted sight lines, components found within 
±300 kms~^ of the targeted galaxy are considered to be as- 
sociated with the galaxy (as described further in i|4]). Those 
components are then grouped into systems. Our definition of 
a system is all the components that lie within 100 kms^^ of 
any other component within that system. We do this because 
fitting detailed line profiles can be non-unique and quite sen- 
sitive to the assumed S/N, as described further in §2.5. The 
temperatures and densities of the systems are calculated as 
the column density-weighted means of the individual com- 
ponents. In most cases, a single component dominates a sys- 
tem. 

The histograms along the temperature and overdensity 
axes indicate the relative fractions of the cosmic mass den- 
sity (grey histograms) and the absorption systems found in 
our lines of sight (red) as a function of temperature and over- 
density. Comparing the red and grey histograms indicates 
how well the absorption seen in each ion at a given impact 
parameter traces the underlying density and temperature 
distribution of all such ions within the volume. Note that 
the histograms are linearly (not logarithmically) scaled, and 
that the integral under the red and grey histograms are set 
to be equal. The grey histograms are the same for all plots of 
a given ion, and sum only gas outside galactic ISM. In con- 



trast, the red histograms vary depending on the absorption 
found for the targeted sight lines at various b. To avoid out- 
liers dominating the computation, the red histograms have 
been generated with a cap in column density of lO^'' cm~^ 
for Hi and 10^^ cm~^ for the metal lines. All the lines above 
these values have been set to this value when summing the 
column density-weighted histogram, to avoid having the his- 
tograms skewed by single large, saturated absorbers, since 
such absorbers generally have highly uncertain column den- 
sities from Voigt profile fitting. This affects 1.3% of all Hi 
absorbers and no more than 0.7% of all metal-line systems. 

The vertical column of panels show the distribution at 
three different targeted impact parameters, 10 kpc, 100 kpc, 
and 1 Mpc, along with random LOS (bottom) for each ionic 
species. The different columns correspond to different ions, 
beginning with Hi and then ordered by increasing ionisa- 
tion potential (discussed below in Figure [S]), namely Mgll, 
Si IV, Civ, Ovi, and Neviii. Each galaxy has four LOS per 
impact parameter, and there are 250 galaxies per mass bin. 
Therefore, each panel of the top three rows represent 1000 
lines of sight. We use a velocity window of ±300 kms~^ for 
the targeted LOS throughout, which we explain in §4 as the 
rough window containing the majority of absorption associ- 
ated with the galaxy. Hence, the total path length in each 
panel is 6 x lO'' kms~^ or Sz = 2.52. In the bottom row. 
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Figure 3. Phase space plots showing the location, in density and temperature space, of Hi (left) and Mgll (right) for = lO'^^ Mq at 
z=0.25. Each coloured point represents an absorption system found along targeted lines of sight at 10 kpc (top row), 100 kpc (second), 
1 Mpc (third), and along random lines of sight (bottom). The grey shading shows the mass- weighted absorption for all gas below 
riH = 0.13 cm~^ outside of all galaxies in the simulation volume, at each location in phase space. This is identical in each plot of a given 
ion. The red and grey histograms show the distributions in identified absorbers and the entire volume, respectively, collapsed along each 
axis. The histograms are linearly scaled, and the integral under the red and grey histograms are set to be equal. Note that the absorber 
colour scale is different for H I. 
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Figure 4. Analogous to Figure \3\ phase spaee plots showing the location, in density and temperature space, of Si IV (left) and C 
(right) absorbers, for Mj, = 10^^ Mq, at the impact parameters indicated. 
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Figure 5. Analogous to Figure [S] phase space plots showing the location, in density and temperature space, of Ovi (left) and Nevill 
(right) absorbers, for = 10^^ Af©! at the impact parameters indicated. 
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we subsample the absorbers from the random lines of sight 
chosen to cover an equivalent redshift pathlength. 

The first clear trend from these figures is that the num- 
ber of systems goes down with increasing impact parame- 
ter. The rate at which the number drops shows some differ- 
ences among the various ions; we will explore this further 
in §6. The straightforward explanation for this is that both 
the metallicity and the gas density drop as one moves away 
from the galaxy, which translates into less metal absorption. 
Nonetheless, it is clear that even at 1 Mpc, there are more 
systems than in the random LOS (see §6). 

Examining the plots more closely, one sees that for some 
ions there is a distinct shift in the overdensities probed by 
that ion as one moves out in impact parameter. This is most 
clear for the higher ionisation potential lines of Ovi and 
Ne VIII, where the peak of the red histograms moves to lower 
overdensities at higher impact parameter, but it is also true 
for all other ions except Mgll. We remind the reader that 
the grey histograms are the same in every panel for each 
species because they correspond to the cosmic mass density 
of that ion, and thus we can interpret the red histograms as 
tracing a subset of the cosmic mass density. 

To illustrate some of the information contained within 
these plots and histograms, consider C iv. At impact param- 
eters of 10 kpc, C IV traces an average nu ~ 10~^ cm~^ (red 
histogram), but the cosmic mass-averaged Civ density is 
riH ~ 10~^ cm~^ (grey histogram). In fact, there is very lit- 
tle cosmic Civ at nn > 10"'^ cm~'^. However, 10 kpc from a 
galaxy inside a 10^^ Mq halo corresponds to a rare location 
in the Universe where there are high levels of metal enrich- 
ment and high densities, and hence Civ traces these metals. 
At impact parameters of 100 kpc, the red histogram for C iv 
peaks at 10"* cm~'^, similar to the cosmic averaged peak in 
the grey histogram. However, the distribution of densities is 
narrower, indicating that absorbers at this impact parame- 
ter do not account for the full range of densities giving rise 
to all the C iv absorption. At impact parameters of 1 Mpc, 
C IV traces overdensities similar to the cosmic average. Fi- 
nally, the random LOS C iv should theoretically trace the 
cosmic distribution of C iv and the red and grey histograms 
should overlap. However, as discussed in lOppenheimer et al.l 
the absorber histogram is biased towards lower den- 
sities. The random LOS sample covering Sz = 2.5 does not 
adequately sample densities with nu > 10"* cm~^ because 
this pathlength is not long enough to probe these rare high 
density regions. This illustrates a central tenet of this paper, 
namely that close-in targeted LOS probe regions of density 
that are not well-sampled via randomly chosen quasar ab- 
sorption sight lines. 

We can similarly examine other metal ions. Random 
LOS showing absorption from Mg ll usually probe gas within 
roughly 10 kpc of galaxies. As a result, random LOS with 
Mgll absorption have the same overdensity as the 10 kpc 
gas. Meanwhile, Mgll absorption from further away gener- 
ally probes lower density gas. For mid-ions, the random LOS 
have significant absorption farther away, and hence 10 kpc 
LOS pick out especially dense gas, but there is little distinc- 
tion between 100 kpc, 1 Mpc, and random. For high ions, 
these arise in more diffuse gas so that the random LOS can 
probe quite low overdensities, and hence the LOS at 10 and 
100 kpc pick out particularly denser gas. 

Finally, for comparison, we examine neutral hydrogen. 



Hi most closely resembles the general behaviours of Mgll, 
in that the cosmic density of H I is heavily weighted toward 
nu > cm~"^ and most cosmic Hi arises at small impact 
parameters going through ISM gas. There are, of course, a 
few key differences between the well-understood behaviour 
of Hi and a metal-line species like Mgll. Most cosmic Hi is 
in damped Lya systems (PL As), gas that is either within 
or near the galactic ISM (| Wolfe et al.ll2005l ). as indicated by 
the steadily rising grey histograms for Hi. The Lya forest 
tracing diffuse gas, i.e. column densities N-^^ ~ 10^^ — 10^^ 
cm~^, corresponds to an insignificant fraction of the total 
cosmic Hi, unlike the high ionisation potential metal lines 
where most of the cosmic absorption arises from IGM and 
CGM densities. Also, Hi can trace a relatively large range of 
phase space because: (i) hydrogen is present at all densities, 
(ii) Lya has a high oscillator strength, and (iii) its pres- 
ence does not rely on metal enrichment. A comparison of 
the red absorber histograms to the grey cosmic density his- 
tograms makes little sense for Hi since all absorbers > lO^'' 
cm^^ have been reduced in our analysis to lO^'' cm~^, which 
hides most of the Lya absorption from damped Lya and 
Lyman-limit systems. However, the red histograms show a 
significant downward shift in overdensity going from impact 
parameter of 10 kpc out to 1 Mpc, reflecting the very differ- 
ent densities probed at these three impact parameters. 

3.3 Physical Conditions vs. Ionisation Potential 

Overall, the phase space plots illustrate a key general trend: 
the higher the ionisation potential of the metal species, the 
lower the overdensity that it traces on average, and hence 
the further away from galaxies it arises. We will argue that 
this is the case in our simulation because these metal species 
that we explore are primarily photo-ionised, except for high 
ionisation potential lines in massive halos. 

To quantify these trends. Figure [H] plots the column 
density- weighted median overdensity (upper panel) - i.e., 
the overdensity above which 50% of the absorption occurs 
- of absorption systems of various ions, as a function of 
their ionisation potential for going from that ion into the 
next higher ionisation state. To compute the column-density 
weighted median value, we take, for each ion, all the systems 
for all the lines of sight and bin them according to their cor- 
responding overdensity (or temperature), as in the red his- 
tograms from figures [31 H] and [S] We then sum the column 
densities of all the systems in a given bin. We plot the over- 
density (or temperature) corresponding to where half of the 
total column density (for all the bins) is above/below that 
value. 

We also plot the column density-weighted median tem- 
perature (lower panel), i.e. the temperature above which 
50% of the absorption occurs. Points are plotted for the 
well-defined median for a sample of 250 galaxies. To illus- 
trate the range, we show vertical range bars (not error bars) 
that span the 16% to 84% enclosing values for the 100 kpc 
bin, as an example. This is the range of the overdensity (or 
temperature) above which 16% of the absorption occurs to 
the overdensity (or temperature) above which 84% of the 
absorption occurs. The different colour points represent the 
three different impact parameters, with black points show- 
ing the global values from the random LOS for comparison. 
These plots are for galaxies in halos with Mhaio ~ 10^^ M©. 
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We plot here the full range of column densities shown in 
Figures [3][5l but these findings do not depend significantly 
on the column density range probed. 

Figure [6] (top panel) shows that the median density of 
metal absorption decreases steadily with ionisation poten- 
tial, but that the rate of the decrease depends on impact 
parameter. Far from galaxies, the decrease is roughly two 
orders of magnitude from Mg ll to Ne vill. This refiects the 
predominantly photo-ionised nature of our absorption lines 
arising in diffuse gas far from galaxies, as elements require 
lower densities to achieve higher ionisation levels. The mean 
overdensity at 1 Mpc is essentially indistinguishable from 
that of the random LOS. At 10 kpc there is nearly as large a 
decrease in the median overdensity as a function of ionisation 
parameter as at 1 Mpc, but in each case the median overden- 
sity is around two orders of magnitude larger at 10 kpc than 
at 1 Mpc. Impact parameters of 100 kpc are an interesting 
counterpoint, showing less than an order of magnitude dif- 
ference in the mean overdensity from Mg ll to Ne vill. As we 
will discuss, this arises because in IO^^Mq halos the high 
ionisation potential lines start to have an additional con- 
tribution from collisionally ionised gas that occurs in the 
denser regions closer to galaxies. 

The horizontal dotted line approximately delineates the 
over density boundary o f virialised halos using Equation 1 
from [Pave et al.l (|2010l ). At 10 kpc, unsurprisingly, all the 
absorption in all the ions arises in gas that is within ha- 
los. At larger impact parameters and for the random LOS, 
the median overdensity depends strongly on the ionisation 
state, with ions having ionisation potentials lower than C iv 
generally arising in halo gas for IO^^Mq halos, while higher 
ionisation potential absorption often arises in gas with over- 
densities lower than that corresponding to galaxy halos. 
Unfortu nately, this absorption is usually weak and photo- 
ionised (jOppenheimer et al.ll2012l '). and hence generally still 
does not trace the so-called missing baryons in the WHIM. 

Turning to the temperatures, Figure [6] shows that the 
median temperature of metal absorption rises steadily with 
ionisation potential. It also shows that, with the exception 
of Nevill, there is essentially no dependence of metal ab- 
sorption gas temperature on the impact parameter. This 
reflects the fact that the majority of metal ions all arise in 
~ ]^q4-4.5 gg^g^ J g photo-ionised gas temperatures, regard- 
less of the impact parameter. Note that the virial tempera- 
ture in this halo mass range is about a million degrees, the 
upper boundary of the plot. 

For Ne vill at large impact parameters, the median tem- 
perature is around 30,000 K, still considerably lower than 
its coUisional ionisation peak temperature of lO'' K, and 
owes to photo-ionised gas. This begins to change as one 
approaches galaxies. At 100 kpc, about half of the absorp- 
tion comes from gas above 200,000 K, and at 10 kpc most 
of the Nevill is hot, with half the absorption coming from 
gas above 500,000 K that is collisionally ionised. Since (as 
is evident from the phase space plots) these Nevill ab- 
sorbers are also stronger, w e obtain a result matching that of 
lOppenheimer et al.l l|2012l l that strong Ne vill absorbers are 
more often collisionally ionised. By using our targeted LOS, 
we show here that this mostly arises from Nevill within 
100 kpc of galaxies, at leas t in th e se 10 ^'^ halos. As hy- 
pothesised in Oppenheimer et al.l (|2012l ). these strong lines 
are probably the only ones that have been studied carefully 



with COS to date ("e.g.. iTripp et al.ll2011^ . but larger and 
deeper samples should uncover a population of photo-ionised 
Nevill arising in more diffuse gas. 

This hot gas close to galaxies is also evident in the Hi, 
where at 100 kpc half the Hi absorption traces gas that is 
above 50,000 K, and much of this gas is above the halo virial 
density. Even at 1 Mpc there is some hot Hi, as the median 
temperature is still fairly high at 40,000 K, even though 
the median overdensity is quite low. As one can see from 
Figure |3l there are a number of absorbers than can arise 
from both hot halo gas and truly diffuse WHIM gas. Hence, 
in principle, broad Hi absorbers could trace truly diffuse 
WHIM gas (i.e. the missing baryons), but one has to be 
careful not to count hot halo gas, which may be detectable 
by other means. 

Our simulations predict that even high ions are mostly 
ph oto-ionised at a ll rad ii. This is at odds with recent work 
bv IStinson et al.1 (|2012l ). whose simulations find that O vi 
is predominantly collisionally ionised in their individual 
galaxy simulations. Their model has the advantage of hav- 
ing higher resolution than ours, but our models have the 
advantage that our feedback prescription has been carefully 
constrained to match a broad range of observations, includ- 
ing IGM enrichment. The main difference is likely that their 
feedback model relies on super-heating gas within the ISM 
that drives hot gaseous outflows, whereas our model follows 
scalings expected for momentum-driven winds in which the 
gas is pushed out via radiation pressure, and therefore is 
not super-heated. Our outflows do heat once they interact 
with surrounding gas, but since they are typically quite en- 
riched, the metal-enhanced cooling rates are rapid. Simply 
put, their feedback model adds hot gas to the halo by con- 
struction, while ours adds cooler gas by construction; it is 
not immediately evident which is closer to correct, and likely 
depends on the details of how winds are actually launched 
which is currently not well understood. IStinson et all (|2012l ) 
did not show the Ovi temperature as a function of radius, 
area-weighted as would be appropriate for an absorption line 
survey, so it is difficult to compare our results directly. But 
these differences highlight that modeling the CGM is not 
a fully solved problem. Examining such detailed statistics 
as line ratios and alignment statistics between O vi and H I 
(and low metal ions, e.g. Sim) provides a way to charac- 
terize the temperature of the O vi gas and discriminate be- 
tween such scenarios. We note that at least in the ran dom 
lines of sight examined in lOppenheimer fc Davj (l2009l ). the 
observed alignment statistics between O vi and H I were bet- 
ter reproduced in a momentum-driven wind scalings model 
as opposed to our constant wind model that yielded more 
colhsionally-ionised O vi. We are conducting such compar- 
isons now against COS-Halos data. 

In summary, different metal ions probe different physi- 
cal conditions around galaxies, with lower ionisation poten- 
tial lines probing denser gas. Mg ll absorbers probe very high 
density gas in and around the ISM of galaxies, while high 
ionisation potential lines probe diffuse gas in the outskirts 
of halos and beyond the virial radius. The absorbing gas 
temperatures generally reflect the photo-ionised nature of 
metal absorption even down to small impact parameters in 
our models, with the notable exception of Ne vill lines that 
can arise in hot gas near galaxies. In the next section we will 
show that O VI in larger halos where a hot CGM is present is 
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Figure 6. Top panel: Median column density-weighted overdensity as a function of ionisation energy for the six ions we consider. Four 
points for each ion are plotted, corresponding to three different b values (all for mass bin 10^2) and random LOS. 100 kpc ( green squares) 
and random (black circles) are plotted at the energy in eV at which the labeled ion is ionised to the next higher level; 10 kpc (red 
triangles) and 1 Mpc (blue diamonds) are shifted for visibility. Points are plotted at the well-defined median; vertical range (not error) 
bars (shown only on the 100 kpc case for clarity) span the 16% to 84% enclosing values. The dotted horizontal line represents the virial 
overdensity at 2 = 0.25. Bottom panel: Similar to the top panel, but for column density-weighted temperatures. 



also mostly collisionally ionised. These results demonstrate 
that spanning a range of ions with low to high ionisation 
potentials can in principle probe a wide range of physical 
conditions in the CGM, but that when a range of ions are 
seen in a single system, it is probably unwise to assume that 
they all arise from the same gas (for the purposes, e.g., of 
CLOUDY modeling). Throughout the rest of this paper, we 
will discuss our results in terms of low (Mg ll) , mid (Si iv and 
Civ), and high ionisation (O vi and Nevill) lines, since this 
provides an underlying physical context for understanding 
the behaviour of these various metal absorbers. 



3.4 Physical Conditions vs. Halo Mass 

The previous section focused entirely on IO^^Mq halos, as 
representative of a typical L* galaxy halo. However, the in- 
creasing coUisional ionisation contribution to the higher ion- 
isation potential species might suggest that there could be 
some halo mass dependence since, as is e viden t in Figure [2] 
and has been shown by e.g. iKeres et al.l l|2005h . large halos 
have substantially more hot gas. 

To investigate any trends with halo mass in Figure[71 we 
plot in the top panels the column density-weighted median 
overdensity and temperature as in Figure |6l here as a func- 



tion of halo mass, focusing on H I (left panel) , O vi (middle 
panel), and Nevill (right panel). As in Figure [S] we plot 
the full range of column densities shown in Figures |3][5l but 
these findings do not depend much on column density range 
probed. We show three different impact parameters, slightly 
offset horizontally for ease of visibility, and here we show the 
16 — 84% range for all cases. The points for galaxies in halos 
of lO^^Mo are identical to those shown in Figure IS] We do 
not show the random LOS here; the values are similar to the 
1 Mpc case. We also do not show the other ions, because in 
those cases there are no discernible trends with halo mass. 

At 10 kpc, for Hi, the absorption arises from low tem- 
perature gas for all mass bins. While Figure [2] shows ha- 
los of mass W^^Mq have high median temperatures. Fig- 
ure [7] shows the HI absorption mostly comes from low tem- 
perature gas present in these halos. For absorption around 
galaxies in halos of 10^^ Mq, most of the Ovi and Nevill 
absorption owes to gas with overdensities lower than that 
found in galaxy halos, since the densities are below the dot- 
ted virial overdensity line. The exception is Ovi at 10 kpc, 
where about half the absorption comes from the outskirts 
of halos. Both high ionisation potential lines have a median 
temperature within the photo-ionised regime at all impact 
parameters (although for Nevill there is a tail to higher 
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temperatures for impact parameters of 1 Mpc). Hence Ovi 
and Nevill absorption around galaxies with halo masses of 
IQ^^Mq actually traces photo-ionised, diffuse IGM gas. This 
occ urs because there is so little hot gas in these small ha- 
los (|Keres et al.l l2005h that there is little opportunity for 
these ions to trace coUisionally ionised gas. 

For galaxies in halos, the temperatures are 

still mostly representative of photo-ionised gas for O vi, but 
Nevill begins to trace somewhat hotter gas, already indicat- 
ing a contribution from coUisionally ionised gas at smaller 
impact parameters. At impact parameters of 10 kpc, this 
gas is within galaxy halos for both ions, and even at 100 
kpc almost all the O vi absorption and about half the Ne vill 
absorption owes to gas with an overdensity consistent with 
being within galaxy halos. At impact parameters of 1 Mpc, 
the gas giving rise to this high ionisation potential absorp- 
tion still mostly arises from gas outside of galaxy halos. 

Finally, for galaxies in IO^^Mq halos, we see a strong 
dependence of physical conditions on impact parameter, 
with the density and temperature both being substantially 
higher at smaller impact parameters. At impact parameters 
of 10 kpc and 100 kpc for both ions, all the absorption owes 
to gas within galaxy halos that is at temperatures indicative 
of collisional ionisation. This reflects the substantial pres- 
ence of hot, virial-temperature gas within these halos. At 
impact parameters of 1 Mpc, however, the absorption still 
mostly arises from photo-ionised, diffuse IGM gas. 

In summary, the high ionisation potential lines tend to 
trace coUisionally ionised hot gas when it is present. But 
this gas is generally only present abundantly well inside of 
halos, at impact parameters of less than 100 kpc, and only in 
massive halos where a hot gaseous atmosphere can form. In 
those cases, these (and other high ion) lines may trace hot 
gas in halos that are not easily probed by X-ray emission 
(or absorption lines), offering a unique opportunity to study 
these baryons. 



4 ABSORPTION AROUND GALAXIES IN 
REDSHIFT SPACE 

The distance of an absorber from a galaxy can be observa- 
tionally characterised by two parameters: the line-of-sight 
velocity difference Aw and the impact parameter 6. In this 
section we discuss the former, i.e. how absorption proper- 
ties vary with the velocity distance from the central galaxy. 
We seek to answer questions like the following: Does the 
absorption drop off with A« at different rates for different 
ions? Can we identify a characteristic LOS velocity distance 
over which the galaxy provides a clear excess of absorption? 
How do these tendencies reflect the physical conditions of 
the absorbing gas? 

Figure [8] shows the column density (N) versus veloc- 
ity separation (An) from the central galaxy for each of our 
six ions, with the species ordered by increasing ionisation 
potential going down the columns. The lines show the me- 
dian column densities for all lines of sight at the three dif- 
ferent impact parameters: 10 kpc (red), 100 kpc (green), 
1 Mpc (blue), for the three different halo masses: 
(dotted), 10^^ M0 (solid), and lO^M© (dashed). We use 
100 kms~^ bins and plot at the bin's midpoint. For example. 
All — 150 kms~^ plots the total column density between 



±(100 — 200) kms~^. Note that the y axis has a rather large 
range; most of this range is not accessible observationally, 
but we include it to accentuate the trends. 

The median column density decreases with increasing 
An for all the ions, following the trend highlighted before 
that the absorption tends to be higher near galaxies. Except 
at high halo masses, there are drops in these curves before 
±300 kms~^. Conservatively, we may say that the vast ma- 
jority of absorption generally occurs within ±300 kms~^ 
of the central galaxy's velocity, although in many cases 
the majority occurs within a smaller velocity interval. We 
choose to define the velocity window associated with the 
galaxy as 300 kms~^ to account for as much absorption 
in the higher mass halos as possible while not going far 
beyond the drops in the lower mass halos. We note that 
±300 k ms~^ is roughly consist ent with recent observational 
work bv lProchaska et al] (|201ll ). Quantitatively, we can con- 
sider the total amount of absorption within ±300 kms~^ 
relative to total absorption within ±600 kms~^. Looking at 
all three impact parameters, O vi has at least 96% of its col- 
umn density within this velocity limit for Mhaio = 10^^ Mq, 
97% for Mhaio = 10^^ M©, and 88% for Mhaio = lO" Mq. 
The absorption of lower ionisation potential metal species 
falls off even faster with increasing velocity difference. This 
indicates that in general one has to look only within the 
central 600 kms~^ window around a galaxy to find most of 
the absorption associated with it. 

Absorption around galaxies in our most massive halo 
bin, IQ^'^Mq, shows a shallower decrease with increasing 
An, reflecting these halos' larger virial velocities and hence 
larger peculiar motions of gas and satellites. While the 
±300 kms"^ cut misses some absorption in the large ha- 
los at small b, LOS passing so close to such massive galax- 
ies are rare. Henceforth, we only consider absorption within 
±300 kms~^ of galaxies for our targeted LOS in the subse- 
quent figures. 

Now we examine the differences between the various 
metal ions. The low ionisation potential ions decrease in 
absorption more rapidly with increasing impact parameter 
than high ionisation potential ions. The column densities 
of Si IV and C iv at fixed An drop rapidly from 10 kpc to 
100 kpc, and at 1 Mpc they are not even visible in this 
plot. In contrast, the high ionisation potential ions, Ovi 
and Nevill, show essentially no decrease in absorption at 
fixed An when the impact parameter increases from 10 kpc 
to 100 kpc, and the drop from 100 kpc to 1 Mpc is more 
modest than for low ionisation absorbers. The decrease in 
column density is also more dramatic with increasing An 
for the low ionisation potential ions. This occurs because 
high ionisation potential lines arise in more extended gas 
distributions at lower overdensities, which trace the general 
large-scale structure in which galaxies live, while the low ion- 
isation potential ions are more confined to the high-density 
gas found closer to galaxies (as seen in Figure [5)). 

In this plot, Hi behaves much like a high ionisation po- 
tential line - note the similarity in the shapes of the curves 
between H l and O Vl, although there is a large difference in 
the magnitude of the column densities. This refiects the fact 
that H I can arise in a wide range of physical conditions, and 
even at small impact parameters there is a substantial con- 
tribution to t he column density from gas along the LOS ou t 
to large An 's jKoUmeier et al.ll2003l : IkoUmeier et aLlbOOd ). 
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Figure 7. Median column density- weiglited overdensity and temperature vs. iialo mass for Hi (left panels), Ovi (middle panels), and 
NeVlll (right panels). For all panels, the different colours correspond to the three different impact parameters: iO kpc (red triangle), iOO 
kpc (green square), and i Mpc (blue diamond) (offset for clarity). The range bars span the 16% to 84% enclosing values, and the dotted 
horizontal line represents the virial overdensity. 




100 200 300 400 500 600 100 200 300 400 500 600 
A V (km/sec) A v (km/sec) 



Figure 8. Median column density (N) vs. velocity separation {Av) from the central galaxy Av for Af^ = IO^^'-'^'^'-'^'^Mq (dotted, solid, 
and dashed lines, respectively). We show results at three impact parameters: 10 kpc (red), 100 kpc (green), and 1 Mpc (blue). Vertical 
dotted lines mark 300 kms~^. 
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At 1 Mpc, Hi absorption does not depend much on the cen- 
tral halo mass. H I is observed almost everywhere that metal 
ions are seen, and in almost every case has a column density 
greater than any metal ion. 

Finally, we examine the trends with halo mass. Nomi- 
nally, the larger peculiar velocities within larger halos should 
result in greater Au's; each factor of 10 in halo mass should 
correspond to a factor of 10^''^ ~ 2.2 in Av. In practise, one 
would expect Aw differences that are somewhat lower than 
this because one integrates the column density through the 
entire halo. Low ionisation potential ions like Mgll are close 
to this expectation, with the differences between the curves 
being roughly a factor of 1.7 in Au. This indicates that Mg ll 
basically arises only when the LOS intercepts a galaxy, and 
at large Ai; the satellite galaxies giving rise to the absorption 
are tracing the underlying halo potential. The mid ionisa- 
tion potential ions Si iv and C iv are similar to Mg ll, in- 
dicating that the denser gas giving rise to these ions also 
traces the underlying potential. Moving towards higher ion- 
isation potential lines, we see somewhat smaller differences 
as a function of halo mass, as absorption in these ions starts 
to pick up gas that is outside the virial radius and hence 
not dominated by the halo potential. It is important to note 
once again that all the massive galaxies in the vzw simula- 
tion have star formation rates well above comparable-mass 
galaxies in the real Universe, owing t o a lack of a quenchin g 
mechanism in our simulations (e.g. iGabor fc Davj [20121 ). 
Therefore, the absorption trends at higher mass may not be 
reflective of the real Universe, where the observed Ovi de- 
clines in strength in more massive halos (|Tumlinson et all 
I2OIII ). Nonetheless, to first order for all ions, close to galax- 
ies it is the dynamics of the host halo that establishes the 
An distribution of the absorption. 

In summary, all ions have the vast majority of their 
absorption arising within a redshift-space distance of (con- 
servatively) Aii± 300 kms~^ around galaxies for all but the 
most massive halos. Low ionisation species show a sharp 
drop with Aw, while high ionisation species show a more 
gradual drop. More massive halos show a broader absorption 
distribution in Aw reflecting their larger potential wells. For 
low ionisation lines, absorption drops very rapidly with im- 
pact parameter, while for high ionisation lines, the drop with 
impact parameter is much slower. These trends broadly re- 
flect the overall morphology of absorption relative to galax- 
ies, in the sense that lower ionisation lines are more conflned 
to dense gas closer to galaxies. 



5 ABSORPTION AROUND GALAXIES 
VERSUS IMPACT PARAMETER 

We now turn to examining the extent of absorption around 
galaxies in physical space, as quantifled by the impact pa- 
rameter b. From the previous section, we know that much of 
the absorption occurs within approximately ±300 km/s of a 
galaxy (the mild exception being for low ionisation potential 
lines in massive halos). Hence, we will examine how absorp- 
tion within this Aw range varies with impact parameter, as 
a function of both ionisation level and halo mass. 

In Figure [5] we plot the summed column density per unit 
redshift along all lines of sight with a given impact parame- 
ter, with the impact parameters ranging from 10— 1000 kpc. 



We show results from our vzw simulation for galaxies in ha- 
los with masses of 1O"M0 (dotted lines), IQ^^A/q (solid), 
and IO'^^Mq (dashed). We show error bars owing to cos- 
mic variance (since statistical errors are small). To calculate 
those, we divide our simulation volume into sixteen sections 
of equal volume and compute the dN/dz values for lines of 
sight within each section, and compute the dispersion over 
the sixteen sections. We ignore the (rare) sections that have 
no absorption when calculating the dispersion. 

Here, we have chosen to plot a summed column density, 
instead of a median one. This is because we are interested 
in predicting the total absorption along the line of sight as 
one moves out in impact parameter. The median value can 
be highly dependent on the resolution and noise level of the 
spectra, since better quality spectra will result in many more 
weak lines. In contrast, the summed absorption is a more ro- 
bust quantity. However, it does have the disadvantage that 
it can be biased by a single, very large absorber. To miti- 
gate this, we apply here the same column density caps as 
described in §3.2. Lines above 10^^ cm~^ for Hi and 10^^ 
cm~^ for metal lines are reset to IQ^^ cm~'^and 10^^ cm~^, 
respectively. While this particular statistic has yet to be de- 
termined observationally as a function of impact parameter, 
this could certainly be done, and would provide a robust 
quantitative estimate of the total amount of absorption as 
one moves away from galaxies. 

The trends seen here in physical space are qualitatively 
similar to those seen in redshift space in the previous sec- 
tion. All the ions show enhanced absorption near galaxies; 
even at 1 Mpc the dNion/dz for targeted LOS is noticeably 
higher than that for random LOS, shown as the diamond 
at 1 Mpc in each plot. For low and mid ionisation potential 
ions, the decline is very steep with increasing radius, and be- 
yond a few hundred kpc, dNion/dz is almost independent of 
impact parameter. For Mgll, the dN/dz value has dropped 
by two orders of magnitude by 100 kpc. For mid ions, the 
sharp drops in dN/dz happen at slightly further impact pa- 
rameters, roughly 200-300 kpc. For high ionisation potential 
ions, the decline is not as steep with impact parameter, but 
there is still a clear enhancement within roughy 300 kpc; 
about a factor of two for O vi and about a factor of four for 
Ne VIII. 

There are also trends with halo mass, although they 
are not strong as in Figure [8] The general trend is that 
there is somewhat more absorption at higher halo masses, 
at most impact parameters. To first order, this reflects the 
increased gas density both in and around larger halos. There 
are some interesting exceptions; for instance, O vi shows 
less absorption within 200 kpc for galaxies in massive ha- 
los. This may reflect temperatures deep within group-sized 
halos that exceed the coUisional ionisation tem perature of 
O VI (jPave. Oppenheimer. fc SivanandamI I2OO8I ) combined 
with densities that exceed the pho to-ionisational densities of 
O VI (lOppenheimer fc Davell2009l ). When one reaches 1 Mpc 
impact parameters, the dNion/dz values for H I and the high 
ionisation potential ions for the various halo masses all have 
essentially converged. 

We have argued above somewhat indirectly that most 
metal absorption arises from gas within roughly 300 kpc 
of galaxies (depending on the ion), but in simulations we 
can test this hypothesis directly. To do so, we tag gas 
particles in our simulations that lie within a chosen ra- 
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Figure 9. The total column density (not number density) per unit redshift versus impact parameter around galaxies in halos of IO^^A/q 
(dotted), 10^^ Mq (solid) and 10^^ Mq (dashed). For halos, dN/dz is the summed column density for the given ion along all lines of sight at 
the given impact parameter over the velocity range ±300 kms~^ , divided by Az = (Number of halos in sample) X (600 km s~^)/(l + z). 
The black symbols indicate dN/dz for random lines of sight. Note that the vertical range for each ion varies, although it always spans 
4 dex. The vertical dotted line corresponds to 300 kpc. Error bars show the cosmic variance across sixteen simulation sections of equal 
volume. 



dius r from any resolved SKID-identified galaxy {M, ^ 
10^'^ Mq), and then regenerate our LOS with contributions 
only from those tagged particles. We exclude galaxies below 
this ma ss limit as they are not well resolved in our simu- 
lations ijFinlator et al.ll2006l ). Fi gure 1101 shows the resulting 
summed dNion/dz values as a function of impact parame- 
ter, analogous to Figure [9] with the different lines showing 
the contribution from particles within r < 30 kpc (blue dot- 
ted), r < 100 kpc (green dashed), and r < 300 kpc (red dot- 
dashed) from galaxies. We show here results for the 10^^ Mq 
halo mass bin, but the results are not significantly different 
for other halo masses. For comparison, the solid black line 
includes absorption from all gas, reproduced from Figure (9] 

Explaining further, in this Figure we use the same 
galaxies as in Figure [9] with the same LOS at the same 
impact parameters. But we no longer consider absorption 
from every SPH particle that intersects the targeted LOS. 
We now consider absorption from the SPIf particles along 
the line of sight only if they are within a given physical ra- 
dius from any galaxy in our simulation volume, not just the 
targeted one. The reader will note that, for example, the red 
300 kpc line extends past 300 kpc. This is due to clustering. 
A targeted line of sight can find absorbers at, for example, 
750 kpc because those absorbers are from particles within 
300 kpc of another galaxy, as opposed to the targeted one. 

For the low and mid metal ions, particles within 300 kpc 
are responsible for the great majority of the absorption at all 
impact parameters, as given in Table (2] For low ionisation 



Table 2. For each ion, we give the percentage of absorption from 
particles within 300 kpc of any galaxy in our simulations, at each 
impact parameter listed. Values in this table are derived from the 
difference between the red dash-dotted and black solid lines in 
Figure fTOl at impact parameters of 10, 100, and 300 kpc. 
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species, particles within even smaller radii are responsible for 
the majority of absorption. For Mg ll, most of the absorption 
owes to particles within only 100 kpc, and typically around 
half of the Mg ll absorption owes to gas within only 30 kpc 
of galaxies. For mid ions Si IV and Civ, essentially all the 
absorption comes from within 300 kpc of galaxies, with the 
majority of it from within 100 kpc, and only a small fraction 
within 30 kpc. For O vi and Ne vill, how much absorption 
comes from particles within 300 kpc depends upon impact 
parameter, as shown in Table [2] For these high ions, very 
little of the absorption comes from particles within 100 kpc, 
and in fact a substantial fraction can come from particles 
further than 300 kpc from galaxies. The contributions to the 
100 kpc curves at 6 > 100 kpc, for example, must come from 
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Figure 10. Total column density per unit redshift (dN/dz) for targeted LOS around IO^^Mq halos, for all gas particles (black solid 
line, identical to Figure |9j, and for gas within the specified radius of a galaxy (coloured lines). Black line is labeled "all", indicating that 
we consider absorption from all SPH particles (i.e. the full simulation) intersecting the targeted LOS. Coloured lines indicate dN/dz 
where we only include absorption from SPH particles that are within a sphere of the indicated physical radius from any galaxy in our 
simulation (not just the targeted one). Dotted purple shows dN/dz for gas within 30 kpc (physical) of a galaxy. Dashed green line is for 
100 kpc, and dash-dot red is for 300 kpc. 



gas associated with other galaxies, either satelhte systems or 
galaxies projected along the LOS with |A«| < 300 kms~^. 

In summary, the extent of metal ions around galaxies 
- that is, the range containing the vast majority of metal 
absorption - is roughly 300 kms~^ in velocity space (al- 
though in some cases much less than that), and roughly 
300 kpc in physical separation (except for Ovi and Nevill 
at large impact parameter). Low ionisation potential metal 
ions are more confined around galaxies in both velocity and 
physical space. This is consistent with the visual impression 
of the column density images in Figure [2] We reiterate that 
metal ions are more likely to be found within 300 kpc or 300 
kms~^ of some galaxy, not necessarily the target galaxy. In 
particular, the excess of low ions at large radii in Figures [9] 
and [To] arises from satellites and neighboring galaxies. 



6 COLUMN DENSITY DISTRIBUTIONS 

The column density distribution (CDD), i.e. the number of 
absorption systems per unit column density per unit red- 
shift, represents the most basic counting statistic for char- 
acterising absorption line systems. Examining the CDDs 
as a function of impact parameter tells us how absorption 
drops off with distance as a function of column density, and 
thereby gives us a more detailed view of how absorption 
varies around galaxies compared to the aggregate statistics 
presented in the previous sections. 



Figure [TT] shows CDDs plotted as N'^ x /(TV) where 
/(TV) = d^n/dNdz for all our ions, where TV is the column 
density, n is the number of lines, and dz is the redshift- 
space path length, corresponding to ±300 kms~^. We plot 
this at three impact parameters of 10 kpc (red), 100 kpc 
(green), and 1 Mpc (blue), and for comparison we also 
plot the CDD for random LOS (black; similar to those in 
lOppcnhoi mer et al.ll2012l ). As before, we combine ions with 
separations Aw < 100 kms~^ into systems before measuring 
the CDD, and apply column density caps (see §3.2). In the 
top panel, we add colored circles indicating the center of the 
distribution, i.e., where 50% of the total column density is 
below that value and 50% above it. In the bottom panel we 
indicate these central values with triangles and squares for 
IQ^^Mq and lO^^TVf©, respectively. The top panels show the 
results for galaxies in halos of ~ IQ^^Mq. The bottom pan- 
els show analogous CDDs for galaxies in halos of ^ lO"Af0 
(dotted lines) and ~ lO^^iVf© (dashed) to assess the depen- 
dence on halo mass. 

We multiply /(TV) by TV^ to obtain a quantity that re- 
flects the (relative) amount of absorption per unit redshift 
at each column density in that ion: the value of TV^ x /(TV) 
at TV = 10^* cm~^, for example, represents the total column 
density per unit redshift contributed by lines in a A In TV = 1 
interval centered at TV = lO^^cm"^. Multiplying by TV^ 
also enhances the visibility of trends by mitigating the typi- 
cally steep power-law dependence of /(TV). In previous work 
lOppenheimer et al ] |2012l : iDave et al.ll2O10l ) , we had plotted 
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Figure 11. Top panel: Column density distributions (CDDs) for galaxies with M^alo = lO^^Af0. /(N) = d^n/dNdz. We multiply /(N) 
by N'^ to show from where most of the absorption per unit log column density arises. The black line is the column density distribution 
for random LOS, red, green, and blue are for targeted LOS at 10 kpc, 100 kpc, and 1 Mpc, respectively. Circles (which increase in size 
with increasing h) mark the peak of the distribution, i.e., where 50% of the column density is below that value and 50% above it. Bottom: 
Same as the top panel for galaxies in halos of IO^^A/q (dotted line) and lO^^AfQ (dashed line). Triangles and squares mark the peak 
of the distribution for IO^^A/q and IO^'^A/q respectively, also increasing in size with increasing b, while circles mark peak for random 
distribution, as in top panel. Vertical lines show the completeness limit of the CDDs for S/N=30 per 6 kms~^ pixel. 
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N X f{N) just for visibility's sake. Here we add an addi- 
tional power of A'', because this means that the crest of the 
A'^^ X f{N) curve represents the column density contribut- 
ing the most absorption per logarithmic interval in that ion. 
These crests are typically shallow — indicating absorption 
that is spread over a fairly wide range of column densities 
— but they correspond well to the centers of the cumulative 
distributions marked by the circles. 

We plot a vertical line that represents the 50% com- 
pleteness limit for the CDD of absorbers identified in our 
artificial spectra with S/N=30. We determine this limit by 
comparing these CDDs to those derived from artificial spec- 
tra with S/N=100 (not shown) and by identifying the col- 
umn density where the S/N=30 CDD begins to deviate by 
more than 50% from the CDD derived with S/N=100. We 
determine this limit using random LOS since this sample 
contains the largest total number of lines owing to its large 
path length. Since we employ the same S/N in all spectra, 
this completeness limit should also be applicable for the tar- 
geted LOS. 

The first thing that one gleans from these figures is that 
for the metals, even at impact parameters of 1 Mpc there 
is clearly more absorption near galaxies than in the random 
LOS. This is true for every ion and for the full range of 
halo masses that we explore in this paper. For low ions that 
arise predominantly near galaxies, this reflects the large- 
scale auto correlation function of galaxies, which extends to 
many megaparsecs. For high ionisation potential ions, the 
absorption arises in less dense (diffuse and WIflM) gas that 
still correlates with galaxies living in large-scale structures. 

We also see that the cosmic absorption in every metal 
ion peaks in the range of « 10^^ - 10^'^cm"^. These 
peaks are generally above the completeness limit for all the 
ions except Nevill, where the lines are intrinsically weak. 
For the lower ionisation potential ions (Mgii, Siiv, Civ), 
this peak moves to larger columns as one goes to smaller 
impact parameters, while for the high ionisation potential 
ions the peak remains mostly independent of the impact 
parameter. In addition, for high ionisation potential ions 
(and for Mg ll) one must probe to fairly low column densities 
to capture the bulk of the cosmic absorption, e.g., A'^ovi ~ 
10^*cm~'^ and A^ncViii ^ lO^^'^cm"^. Hence, while many 
O VI and Ne v ill absorption systems are now being det ected 
in COS data l|Tumlinson et al.ll201ll : [Tripp et al.ll2011^ ■ our 
models predict that higher S/N observations that probe to 
A'^ ~ lO^^'^cm"^ are needed to capture the bulk of cosmic 
absorption in these ions. 

Now let us examine Hi, as this ion is relatively well un- 
derstood in terms of i ts connection with galaxies and large - 
scale structure (e.g., iDave et al.l Il999l : iDave et all I2OI0I ). 
H I shows a clear trend of having an increasing number of 
high column density absorption systems as one approaches a 
galaxy, i.e. as one goes to small impact parameters. In con- 
trast, as one goes to smaller impact parameters the number 
of lower column density absorption systems does not change 
as significantly. In fact, below N ~ lO^'^cm"^ the trend ac- 
tually reverses with the number of systems increasing as 
one goes to larger impact parameters. This can be easily 
understood. The strong absorbers mostly arise in the dense 
gas around galaxies, while the weak absorbers arise more 
in the surrounding large-scale structure, even when the im- 
pact parameter is small. This is consistent with the strong 



correlation between H I colu mn density and ov erdensity pre- 



dieted in the models (e.g., Dave et al] I2001I : [Schava i20()ll : 



iDave et al.ll2010l : iHui et al.Hl997l ) 

For the metal lines, the variation in the shape of the 
CDD with impact parameter shows interesting trends as a 
function of ionisation level. As one moves from impact pa- 
rameters of 100 kpc to 10 kpc, for lower column density 
absorbers there is little change in the incidence of absorp- 
tion. However, as one moves to higher column densities, the 
number of absorbers increases and the peak moves to higher 
column densities. The column density where the increase 
to small impact parameter becomes noticeable increases 
with ionisation level. For instance, for Mgll the 10 kpc and 
100 kpc curves differ at A^Mgii > lO^^cm"^, while for Civ 
they differ at A'civ > lO^^cm"^. For Ovi and Nevill, there 
are only minimal differences in the CDDs at impact param- 
eters of 10 kpc and 100 kpc. At impact parameters of 1 Mpc 
the number of absorbers is smaller at all columns for all the 
metal ions, with the peak also generally occurring at lower 
column densities. The random LOS continue this trend. In 
short, the lower the ionisation potential, the more sensitive 
the ion's CDD is to the proximity of a galaxy. 

It is conventional to fit the CDD with a power law in col- 
umn density for a given ion: f{N) oc . For random lines 
of sight, the observed power-law slopes for H I, O VI, and C IV 
are generally b etweei i 1.5 < /3 < 2.2 jDanforth et al.|[200^ : 
ICooksev et al]|2010l '). lDave et al]||2010i r using similar simu- 
lations, found a slope for the Hi CDD of /3 = 1.70, which is 
(unsurprisingly) similar to what we find here. The slope for 
O VI depends on the range of column densities over which the 
fit is done, since a pure power law is not a good descriptor. 
Using the full range above our c ompleteness limi t , we o btain 
/3 ~ 2.3, which is comparable to lDanforth et al] (|2006l ) who 
found P — 2.2±0.1. Civ likewise is not a perfect power law, 
but we find a slope of /3 ~ 1.9, w hich is steeper that ob- 
servations bv lCooksev et al.l (|2010l ) that yield ^ = 1-5^0.11, 
but probably within uncertainties given that incompleteness 
in the data which has poorer quality than our simulated 
spectra will generically lead to shallower slopes. Given the 
variations in spectral quality between all these data sets, we 
consider our simulations to be broadly in agreement with 
current measures of CDD slopes for these ions. Observations 
have yet to constrain the slope as a function of impact pa- 
rameter, but our models predict that the slope does not vary 
dramatically, and mainly only the amplitude increases sub- 
stantially to small impact parameter. This prediction will 
be testable with upcoming observations that have sufficient 
statistics to examine f{N) as a function of impact parame- 
ter. 

Looking at the bottom panels of Figure 1111 absorption 
around galaxies in halos of larger and smaller masses con- 
tinues these same general trends. At impact parameters of 
10 kpc, as one goes from halo with masses of IO^^AIq to 
10^'^ A/0 one sees similar trends as when one went from im- 
pact parameters of 100 kpc to 10 kpc in halos of fixed mass; 
there is agreement of the CDDs at small column densities 
and more absorption at high column densities, and the tran- 
sition occurs at a column density that increases as the ion- 
isation potential of the ion increases. Hence in this regard, 
going to higher halo masses is equivalent to going to smaller 
impact parameters in a halo of fixed mass. The CDDs are 
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relatively insensitive to the mass of the galaxy halo for im- 
pact parameters of 1 Mpc. 

Ne Vlll exhibits an interesting trend at high halo masses, 
arising from its strong coUisional ionisation contribution 
when a hot gaseous atmosphere is present. We see in Fig- 
ure [S] that Ne vill has a low-density, photo-ionised compo- 
nent probing the diffuse IGM, which gives it the more ex- 
tended characteristics of a high ionisation potential ion like 
O VI that is itself mostly photo- ionised. However, Ne vill also 
traces iq^-^-^ o k hot halo gas, which is both denser and at 
lower impact parameters, where cooler ions like Siiv and 
C IV are also found. Therefore for large halos and impact 
parameters inside of 100 kpc, there is a substantial popu- 
lation of (coUisionally ionised) Nevill absorbers. Note that 
this Ne VIII could be coincident wit h lower ionis a tion p oten- 
tial ions, as has been observed by iTripp et al.l (120111 ). but 
arises in a different gas component. 

In summary, the ODD of If I and metal ions shows 
trends with impact parameter that reflect correlations of ab- 
sorption with both nearby galaxies and large-scale structure. 
All ions show increased absorption closer to galaxies. Low 
ionisation potential ions are more influenced by the presence 
of nearby galaxies. These CDDs represent a prediction of hi- 
erarchical models that enrich the IGM using outflows from 
star-forming galaxies and can i n principle be tested an d con- 
strained by observations (e.g.. pTumlinson et al.|[2"oll^ . The 
full COS-Halos data set is being analyzed now (Werk et al., 
2012), and in future work we will undertake a detailed com- 
parison to CDD and other observed statistics, including an 
artiflcial spectra sample that more closely mimics the COS- 
Halos spectra. 



7 VARIATIONS WITH OUTFLOW MODEL 

One expects the enrichment of the surrounding CGM gas 
to depend sensitively on the properties of the enriching out- 
flows. So far we have only considered our favoured outflow 
model with momentum-driven wind scalings. Here we con- 
sider how the absorption properties around galaxies depend 
on our assumed outflow model, using two other wind pre- 
scriptions: a simulation with no winds (nw) and a constant 
wind (cw) model where we assume a constant mass loading 
factor of r; = 2 and a constant wind speed from all galaxies 
of Vw = 680 kms~^. The latter is similar to that used in 
the Over whelmingly Large S imulations (OWLS) reference 
model of ISchave etal] (|2010l '). The cw version used in this 
paper is identical to that des cribed in iDave et al.N2010f) . ex- 
cept that it was re-run using |w iersma et al.l (|2009al ) met al- 
line cooling rates, for consistency with the vzw model used 
here. 

To begin, Figure [2] shows a pictorial representation of 
how the physical conditions and absorption vary with wind 
model. In Figure [T^ for easier comparison, we collapse the 
density and temperature information given in Figure[2]down 
into one dimension by taking the azimuthal average of those 
images. In Figure \T2l we show a larger region than in Fig- 
ure[2]to illustrate larger-scale trends, and include metallicity 
to more fully understand the cooling processes. 

In the top panel of Figure [T^ vzw curiously shows more 
similarities to the no wind model (nw) than to the constant 
wind model (cw). This illustrates that the high wind speeds 



from all galaxies in the cw model cause signiflcantly more 
spatial dispersal of mass (along with metals, as we show be- 
low) on ~Mpc scales around galaxies. The vzw model has 
lower velocity winds that do not have such a dramatic im- 
pact on large scales, but do have a strong impact on smaller 
(CGM) scales. The vzw model pushes more mass via winds 
into the CGM relative to the no wind model, but this mass 
is close enough that it re cycles back onto the galax y in much 
less than a Hubble time (|Oppenheimer et al.ll2O10l ). We note 
that the galaxy population is broadly more similar in the 
two wind models, as they both suppress global star forma- 
tion substantially relative to the no-wind case (for a fuller 
discussion of thes e properties, see lOppenheimer et al.|[201ol : 
iDave et al.ll2011al lbli. but Figure [T2l shows that they are quite 
different in where they deposit the ejected material. 

The middle and lower panels of Figure [12] illustrate how 
winds affect the CGM. Here vzw and cw actually trend in 
the opposite direction relative to the no-wind case! In the 
vzw case, the area around the galaxy is slightly colder than 
with no winds. This is because the lower wind speeds deposit 
more metals around galaxies (lower pa nel), and this results 
in an increased amount of metal cooling l|Oppenheimer et al.l 
that more than offsets the shock heating from the 
winds. In contrast, the cw model expels gas at high veloc- 
ities, around the escape velocity for 10^^ A/0 halos, even 
from small galaxies. This means the enrichment is more 
widespread - past 100 kpc, cw shows greater metallicity 
than the vzw or nw model. Moreover, the wind energy is de- 
posited into less dense gas where it does not have a chance 
to radiate away its energy. 

The H I maps shown in Figure [2]illustrate that winds do 
not make a large difference to H l absorption, at least on the 
large scales depicted here, similar to results at hi gh redshift 
(e.g. , iKoUmeier et al.ll2003l : iKoUmeier et al ] |2006l ). However, 
there is an increase in Hi on small scales, roughly 100 kpc. 
The metal absorption, in contrast, shows more dramatic dif- 
ferences between the wind models. In the no-wind case, all 
the metals are essentially conflned to be in and around galax- 
ies, showing that even the distribution of metals within the 
CGM requires winds. These differences should be manifest 
in the statistics of absorbers around galaxies, providing an 
opportunity to constrain wind models. 

To further quantify the extent of the metal distribution, 
Figure [13] shows the fraction of all cosmic metal mass that 
lies within a given radius from galaxies with Mt ^ 10^'^ Mq 
in the simulation volume. We also plot the fraction of all 
cosmic mass for all species, in addition to just metals. The 
procedure for this was discussed in [JS] Here the impact of 
winds on distributing cosmic metals is shown clearly: the no 
wind model keeps essentially all metals confined very close 
to galaxies, while the constant wind model disperses them 
over large scales, with the momentum-driven wind model 
intermediate between the two. This figure also shows that 
the constant wind model also pushes the total mass further 
out, while the momentum-driven wind and no wind models 
are more similar. For our favoured momentum-driven wind 
scaling model, 83% of the metal mass is in gas within 300 kpc 
from galaxies. For the cw model, only 40% is within this 
radius. Hence the basic extent of metals around galaxies 
can already provide a discriminant between wind models. 

We now aim to quantify these differences in metal dis- 
tribution using absorption line statistics. Figure [14] is similar 
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Figure 12. Median overdensity, temperature, and metallicity vs. impact parameter for three different wind models; momentum-driven 
"vzw" , constant wind "cw" , and no wind "nw" . Sight lines are all around galaxies with halo masses of 10^^ Mq. The top two rows are 
similar to the first two columns of Figure [2] but collapsing the image down to one dimension. 
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Figure 13. Thin lines: Metal mass fraction of all gas particles 
within spheres of radius r around all galaxies Af, > 10®'^ Mq in 
our simulation at z=0.25, relative to the total metal mass of gas 
particles in the whole simulation volume, for momentum-driven 
winds (solid line), no winds (dashed line), and constant winds 
(dotted line). The vertical dotted line delineates 300 kpc. Thick 
lines: Total mass fraction for all species, not just metals. 



to Figure [9] except that it compares the three wind models, 
using the same column density caps as explained in §3.2. 
Here, we also include a comparison to our old vzw simula- 
tion that used collisionally ionised equilibrium (CIE) metal 
line cooling, labeled as vzw-cie, as opposed to our current 



vzw model that uses thelWiersm a et all (j2009ah (PIE) cool- 
ing. For the most part, the differences between the vzw-cie 
model and our current vzw model are small compared to the 
differences between wind models, and hence this aspect of 
our modeling does not introduce a large uncertainty into the 
results. The differences between vzw and vzw-cie are small 
because the winds have moderate velocities, preventing large 
amounts of gas to reach lower overdensities. Instead, they re- 
main in moderate-density gas where cooling times are short 
regardless of photo-ionisation suppression. Note that when 
we generate spectra, we compute ionic abundances including 
photoionisation in both cases; it is only during the evolution 
of the simulation that vzw-cie is different. 

In this figure, as in Figure |5] we also show dNion/dz 
for random LOS for each wind model as the diamonds near 
the right edge. For all models and ions, dNion/dz values at 
1 Mpc are still higher than dNion/dz for random LOS. Note 
that for all the metal ions, there is so little absorption in 
the no wind model that the random LOS dNion/dz (purple 
diamonds) falls off the bottom of this plot. Similarly, the cw 
model shows so little Siiv and Mgll absorption in random 
LOS that the red diamonds fall off the bottom of the plot. 

For all the metal ions, the no wind model gives signif- 
icantly less absorption than any of the wind models, even 
down to the smallest impact parameters probed here. How- 
ever, there is not a large difference between models with 
winds and the model without winds for Hi. Without out- 
flows, metals basically exist outside the ISM only owing 
to tidal or ram pressure stripping processes that remove 
material from the ISM of the central or satellite galaxies. 
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It is evident from this plot (as with the images in Fig- 
ure I12p that such stripping processes provide only a small 
contribution to the CGM metal absorption in these halos, 
though it can be more substantial in large halos (e.g . , 
iDave. Qppenheimer. fc Sivanandanil l2008l : IZu et all I2OIII ). 
Hence, we predict metals seen at any impact parameter be- 
yond that of the ISM of typical galaxies arise almost exclu- 
sively owing to outflows. 

The differences between the constant wind and 
momentum-driven wind cases are more subtle. The global 
trends with ionisation state are as before: lower ionisation 
potential species are more highly peaked at small impact pa- 
rameters in both wind models. The main significant differ- 
ence is an overall offset. In general, the cw model has lower 
absorption for low ionisation potential species, and higher 
absorption for high ionisation potential species relative to 
the vzw (or vzw-cie) model. The lower absorption arises in 
part because the cw model produces les s metals overall ow- 
ing t o its reduced cosmic star formation l|Oppenheimer et al.l 
I2OI2, '). Going to higher ionisation species, absorption in the 
cw model becomes only slightly less than absorption in the 
vzw model. For Nevill, cw shows more absorption than vzw. 
This reflects the impression from Figure [12] that the cw 
model expels metals to greater distances and heats this dif- 
fuse gas more l|Qppenheimer et al.ll2012l ). Hence, the global 
absorption strength as a function of ionisation state provides 
another potential discriminant between outflow models. 

In Figure 1151 we plot the median column density ver- 
sus velocity, analogous to Figure [8] except now we vary the 
wind model as opposed to the halo mass. Not surprisingly, 
for Hi all the wind models give similar results; however vzw 
does show an increase over nw and cw at 100 kpc, because 
these winds push more cool gas to CGM distances. For the 
metal lines, once again the no-wind case shows virtually no 
absorption except perhaps very (dynamically) close to galax- 
ies (±100 kms~^). Hence, tidal effects and other stripping 
processes are ineffective in distributing metals in velocity 
space, just as in physical space. 

In this figure, we see similar trends with wind model 
as we did in Figure 1141 For low ionisation potential metal 
species, the vzw model produces higher column densities, 
while the opposite is true for high ionisation potential 
species. Previous work (jOppenheimer et al. 2012) has shown 
that the cw and vzw models enrich the dense gas in quite 
different ways. The cw model, with its high velocities ema- 
nating even from small galaxies, deposits fewer metals into 
the high density regions very close to galaxies and more met- 
als into the diffuse IGM. Accordingly, metal ions that show 
more absorption in the cw model than in the vzw model 
have more absorption from more diffuse IGM gas. 

In summary, the rate of the decrease in absorption 
of metal lines as one moves away from the central galaxy 
provides a potentially strong discriminant between outflow 
models that enrich the diffuse IGM to fairly similar levels. 
Winds are required to enrich the CGM of normal galaxies 
to any significant level. Our currently favoured momentum- 
driven wind scaling model predicts more low ionisation po- 
tential absorption close to galaxies than the constant wind 
model and less high ionisation potential absorption farther 
from galaxies. Quantitative comparisons with present and 
upcoming COS data should yield more stringent constraints 
on outflow propagation. 



8 CONCLUSIONS 

We have examined the absorption line properties of H I and 
flve key metal ions in the vicinity of galaxies at z — 0.25 in 
cosmological hydrodynamic simulations that include galac- 
tic outflows, as a function of velocity separation, impact pa- 
rameter between the galaxy and the line of sight, and halo 
mass. Our chosen metal ions span a range of ionisation po- 
tentials, from low (Mgii) to mid (Siiv & Civ) to high (Ovi 
& Nevill). This is the first absorption line study of metal 
absorption around galaxies in cosmological hydrodynamic 
simulations with a self-consistently generated outflow model 
that matches a wide range of galaxy and IGM observables. 
Our work is motivated by the upcoming wealth of data from 
Hubble 's Cosmic Origins Spectrograph, and hence our artifi- 
cial spectra are generated with resolutions and noise charac- 
teristics comparable to the best data that will be obtained 
with COS. 

Our primary conclusions are: 

(i) Absorption in all ions is enhanced closer to galaxies. 
This is true in terms of both impact parameter and veloc- 
ity separation. A velocity separation of ±300 kms~^ around 
galaxies encompasses most of the cosmic metal line absorp- 
tion (although in some cases it can be less), and absorption 
is also significantly stronger within 300 kpc around a galaxy, 
particularly for lower ionisation potential lines. These ranges 
are somewhat larger around bigger galaxies. In our favoured 
momentum-driven wind simulation, ~ 80% of all cosmic 
metals lie within 300 kpc of a galaxy. 

(ii) The dependence of metal absorption strength on dis- 
tance from a galaxy in either physical or redshift space 
depends monotonically on the ionisation potential of the 
absorbing ion. Ions with a low ionisation potential (Mgll, 
Silv) arise in higher density gas that tends to drop off 
more quickly with impact parameter and velocity separa- 
tion than ions with higher ionisation potentials (C IV, O vi, 
Nevill). High ionisation potential absorbers are more asso- 
ciated with gas at halo-like overdensities, and the overall 
large-scale structure that contains galaxies, while low ion- 
isation potential lines arise more in the dense gas close to 
individual galaxies. Regardless of ionisation level, even out 
to 1 Mpc, targeted LOS show an excess of absorption over 
random LOS, reflective of large-scale matter clustering. 

(iii) The majority of cosmic absorption in the ions con- 
sidered here occurs in photo-ionised gas at T < 10^ K. The 
exception to this is for the high ionisation lines O vi and 
Nevill at impact parameters < 100 kpc arising in massive 
halos containing substantial hot gas; these lines are predom- 
inantly coUisionally ionised. 

(iv) The dependence of the column density distributions 
on impact parameter also shows trends with ionisation po- 
tential. Species with lower ionisation potentials are more 
affected by the proximity of a galaxy. Excepting Ne vill, our 
spectra (comparable to COS data quality) directly trace the 
majority of cosmic absorption; i.e. there is not a large popu- 
lation of smaller absorbers that is inaccessible to COS that 
would dominate the total cosmic absorption. 

(v) Without winds, even gas within « 100 kpc of galax- 
ies remains mostly unenriched; hence outflows are required 
to enrich the CGM as well as the IGM. The differences in 
metal absorption between our favoured momentum-driven 
wind scaling model and a constant wind model are evident 
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Figure 14. The total column density per unit rcdshift for the vzw, vzw-cie, cw, and nw models, all for galaxies in IO^^Mq halos. The 
diamonds indicate the dNion/dz value for random lines of sight for each wind model. The symbols are shifted slightly to the left of 1 
Mpc, and separated for easier viewing. In some cases the symbols are lower than the scale of the plot so do not appear. 
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Figure 15. The median column density (N) vs. velocity separation from the central galaxy (Ai)) for galaxies in lO^^M© halos for the 
vzw, cw, and nw models. We show results at three impact parameters: 10 kpc (red), 100 kpc (green), and 1 Mpc (blue). The vertical 
dotted marks 300 kms~^. 
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both for lower ionisation potential species close to galaxies 
and for high ionisation potential species farther from galax- 
ies. The models are clearly distinguishable from each other 
using the combined CGM statistics of O vi and Hi. 

These results provide a starting point for understand- 
ing how absorption lines trace the metal enrichment around 
galaxies. We have highlighted some basic trends, and shown 
that quantifying this distribution as a function of impact 
parameter and velocity separation can provide interesting 
constraints on key physical processes such as galactic out- 
flows, as well as the density and temperature state of the 
metal-enriched IGM. 

While we have considered only five metal ions in this pa- 
per, the basically monotonic trends with ionisation potential 
suggest that the behaviour of any other ion can be predicted 
just based on its ionisation potential. Our simulations pre- 
dict that higher ionisation potential ions should fall off more 
slowly at larger impact parameters and velocity separations, 
because they are tracing lower density, mostly photo-ionised 
gas. This suggests a physical structure of CGM gas in which 
O VI is more extended than C iv, which is more extended 
than Si iv, and so on. Ne vill and O vi, if probed at suf- 
ficiently low column densities, trace the most remote and 
diffuse metals of any UV resonance line, but their strongest 
lines arise in hot collisionally ionised gas near large galaxies. 

This work is the first step in a series of works to confront 
successful models for galaxy-IGM coevolution with absorp- 
tion line observations around galaxies in the low-z universe. 
We examine redshift z = 0.25, but the basic trends are ap- 
plicable to all redsh ifts probed by COS (z < 1), as recent 
work bv iDave et al.i t201Q ) and lOppenheimer et al.l (|2012l ') 
show little evolution in the IGM from z = 1 — > 0. We cau- 
tion, though, that detailed comparisons to observations are 
premature for the simulated spectra presented here, since 
quantitative trends can be sensitive to the details of spectral 
resolution, noise level, etc. Furthermore, our massive galax- 
ies (Mhaio ~ 10^^ Mq) are all star- forming in these simu- 
lations, in clear conflict with observations of mostly passive 
galax;ies in this halo mass range, and hence the predictions 
for such halos may be influenced by physical effects not in- 
cluded in our current models. Nonetheless, we believe that 
the basic intuition of understanding metal absorption sur- 
rounding galaxies in terms of the ionisation potential of the 
tracer ion is robust, and provides a clear intuition for in- 
terpreting current and future observations. As observations 
progress, particularly with COS, simulations like these will 
provide a critical testbed for galaxy formation models, and 
will help elucidate the physical processes that drive the en- 
richment of the intergalactic and circumgalactic medium. 
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